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ABSTRACT 

A line list of vibration-rotation transitions for ^'^C substituted HCN is presented. The 
line list is constructed using known experimental levels where available, calculated 
levels and ah initio line intensities originally calculated for the major isotopologue. 
Synthetic spectra are generated and compared with observations for cool carbon star 
WZ Cas. It is suggested that high resolution HCN spectra recorded near 14 ^m should 
be particularly sensitive to the ^'^C to ^^C ratio. 
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1 INTRODUCTION. 

Carbon giant stars are though to arise from the third dredge 
up of asymptotic giant branch (AGB) stars, which pollutes 
the envelope and atmosphere with nuclear processed ma- 
terial from the interior. This increases the abundance of 
carbon and the ^^C/^^C ratio to well above the terrestrial 
level (0.011). In fact the ^''C/^^C ratio in car bon giants has 
been measured to be as high as one-third (|Abia fc IsernI 



1 19971 ). There are a number line lists or opacity functions 
available for H^^ CN (|j0rgensen et al.lll985l : |j0rgensenlll99Ol : 
lAoki et al.lll998l : [Harris et al.ll2002bl . l2006l '). However, when 
considering the contribution to opacity by molecular species 
such as HCN it may not be sufficient to account only for 
H^^CN, but also for the isotopologue H^^CN. 

The calculation of a complete triatomic linelist is com- 
putationally exp ensive, req uiring several tens of thousands 
of CPU hours ( jTennvson e t al. 2007). However, within the 
Born-Oppenheimer approximation the electronic structure 
of isotopologues is identical. Thus both the potential en- 
ergy and electric dipole moment functions are identical for 
all isotopologues. This implies that the vibration-rotation 
frequencies and transition intensities are likely to be very 
similar. 

For a hetronuclear molecule, such as HCN, the main 
differences between the spectra of isotopologues are caused 
by the change in reduced mass. In the harmonic approxi- 
mation the vibrational contribution to the line frequency is 
proportional to p~^ , where p — mim2/{m\ +7712) is the re- 
duced mass. So the frequency of the CN stretch vibrational 



mode of H^^CN and H"CN will differ by the order of 1% 
and the bend and H-C stretch mode by less. In terms of 
both the observation of H^^CN and the line blanketing of 
a model atmosphere this shift in line frequency is more im- 
portant than the small differences in line intensity between 
comparable lines of H^^CN and H^^^CN. 

It is the aim of this work to compute a set of energy 
levels and line frequencies for H^'^CN and HN^^C. These en- 
ergy levels will be used in conjunction with th e Einstein A 
coefficients for H^^CN and HN^'^C computed bv lHarris et al.l 
(2002^), to generate a H"CN/HN^^C linelist. The line fre- 
quencies in this H^'^CN/HN^'^C linelist have been corrected 
using available laboratory data in order to increase its accu- 
racy. The new linelist has been used to compute synthetic 
spectra for C-stars, with different ^'^C/^'^C ratios using star 
WZ Cas as a prototype. 



2 CONSTRUCTION OF A H^'CN AND HN^'^C 
LINELIST. 

2.1 ab initio energy levels. 

Using the existing HCN /HNC potential energy surface of 
Ivan Mourik et al.l (|200in a set of ab initio rotation- vibration 
energy levels for H^'^CN and HN^'^C were calculated for an- 
gular momenta of J = 0, 1, 2, 3, 5, 10, 20, 30, 40, 60 and for 
both even (e) and odd (/) parity. The states were computed 
up to an energy of at least 10000 + -B(J( J + 1)) cm"'^ above 
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the H^^CN ground state, where B ~ 1.5 cm~^ is the HCN 

rotational constant. 

The DVR3D suite of codes IX ennvson et al.ll2004 ) was 
used for all these calculations. DVR3D uses an exact kinetic 
energy operator and discrete variable representation for the 
nuclear motion wavefunctions. We hav e used the init i al basis 
set functions and parameters used bv lHarris et al.l l|2002bl ) 
for the computation of an extensive H^^CN/HN^'^C linelist. 
The basis functions are Legendre polynomials for the angu- 
lar grid points and Morse oscillator-like functions for the ra- 
dial grids. Jacobi coordinates were used, with 50 grid points 
for the angular coordinate, 35 grid points for the first (_R) 
radial coordinate, 21 for the second (r) radial coordinate. 
Where r is the distance from C to N nuclei and R the dis- 
tance from the H nucleus to the centre of mass of the C-N 
diatom. The parameters for the Morse oscillator like basis 
functions in the r coordinate are r^ ~ 2.3 ao, De{r) = 29.0 
Eh and uje{r) = 0.0105 Eh and Re = 3.2 ao, De{R) = 5.0 
Eh and uJe{R) = 0.004 Eh for the R coordinate. Where re is 
the equilibrium distance, De is the dissociation energy and 
uie is the harmonic frequency. 



2.2 Vibrational assignments and rotational 
constants. 

Only the exact quantum numbers for a hetronuclear tri- 
atomic, angular momentum (J) and parity are known for 
each state computed by DVR3D. To determine the rota- 
tional constants for each vibrational state it is first necessary 
to assign approximate vibrational quantum numbers to each 
rotation- vibration state. The approximate vibrational quan- 
tum numbers are C-H stretch (ui), bend (^2), C-N stretch 
(ua), vibrational angular momentum (l) and isomer (H'^^CN 
or HN^^C). Where I take values in steps of 2 from or 1, if 
V2 is even or odd respectively, up to the lower of V2 or J. So 
for a state with 112 = 5 and J ^ 5 then I can take values of 
1,3,5, but if J = 3 or J = 4 then I can only take values of 1 
or 3. 

The quantum number assign ments were made us ing a 
method similar to that described in [Harris et al.l (|2006l ). Ini- 
tially quantum numbers were assigned to the H^^CN J = 
purely vibrational states up to 10 000 cm~^ above the zero 
point energy. The first states to be assigned were the lowest 
lying states of each of the 3 H'^^CN fundamental modes. The 
vibrational expansion eq.[T]was then least squares fit to the 
assigned states, which provided values for the fundamental 
frequencies Ui. 

3 
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The fundamental frequencies allow estimates for the en- 
ergy of more highly excited states to be made, which aid 
further assignments. Next, eq.[T]was least squares fit to the 
newly assigned states so that the higher order vibrational 
constants Xij could also be determined. This process was 
repeated until aU J = H"CN states up to 10 000 cm"^ 



were assigned. The J=0 HN^^C states were fit separately 
from the H^^CN states, but the assignments made simulta- 
neously. The J — 1 e and/ states were assigned by compari- 
son with the J = states and by fitting equation[T]to verify 
the assignments. 

In order to aid the assignments of states with J > 1 we 
have least squares fit the assigned levels for each vibrational 
state and each parity with: 

E{v,J) = E{v) + B4J(J + l) + l^] 

-D4J{J+l) + l'f + ... (2) 

where E{v) is the vibrational energy of the band, Bv is the 
rotational constant and Dy is the centrifugal distortion con- 
stant. This allows energies for the states with the next high- 
est J to be estimated and assigned. For further verification 
these assigned states were then fitted with the vibrational 
expansion. This process was carried out on an increasing J 
by J basis until all the computed energy levels up to J = 60 
had been assigned. Once the assignments were completed, a 
final fit of the rotational expansion was performed for each 
vibrational state and parity. A total of 390 vibrational states 
of H^^CN/HN^^C have been studied all with a vibrational 
energy {E{v)) less than 10 000 cm~^. Using this final set 
of rotational constants a set of calculated rotation- vibration 
energy levels for H^^CN and HN^'^C was computed for all 
values of J up to 60. 

At stellar temperatures transitions between levels with 
lower state energy over 10 000 cm~^ do contribute signifi- 
cantly to opacity. These lines tend to be weak and numerous 
so that they form an almost continuous opacity source. In 
order to account for these high temperature weak lines we 
have supplemented the H^^CN/HN^'^C energy level list with 
high energy H^^CN/HN^^C energy levels. 



2.3 Energy levels determined from laboratory 
data. 

Line frequencies determined in the laboratory are signifi- 
cantly more accurate than than purely ab imtio data. There- 
fore, where possible we have incorporated laboratory deter- 
mined energy levels into our linelist. H'^^ CN has be en well 
studied in the laborator y ( Lehmann et al.lll98 2: Smi th et al.l 
19891: iMaki et all |l_995; Ma ki et al.l bOOotlPevi et al.l 120031 



2004), but HN^^C has been studied to a far lesser extent 



( Maki fc Mellaull2001^ . HN"C laboratory determined en- 
ergy levels have therefore not been incorporated in to the 

linelist. 

iMaki et al.l (|2000l ) have provided an electronic database 
of H^'^CN line frequencies, from this we have compiled a list 
of H^'^C N laboratory determined energy levels, in the same 
way as [Harris et al.l (120061 ). Many line frequencies of the 
studied bands remain unmeasured, there are therefore miss- 
ing energy levels for many of the vibrational states. These 
missing energy levels were interpolated by means of the ro- 
tational expansion to give a complete list list of energy lev- 
els in J for each band up to the maximum J measured by 
iMaki et al.l (|2000l ) . As the rotational expansion is divergent, 
it cannot be reliably used to extrapolate the laboratory de- 
ter mined ener g y lev els beyond the maximum J measured 
bv lMaki et al.l (2000|). To extend the laboratory determined 
energy level list to J = 60, the maximum extent of the 
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[Harris et alj l|2002bl ) linelist, we have used a correction to 
the ah inito energy levels described in 12.21 This constant is 

given by: C = Elab{Jmax) — Eai{Jmax), where Elab{Jmax) 

is the energy of the state in the laboratory energy level list 
with the highest angular momentum {Jmax), and Eai(Jmax), 
is the energy of the ab initio energy level with J = Jmax- 
This lab/empirical energy level list contains 4425 energy lev- 
els, which have been incorporated into the calculated list of 
H"CN/HN"C energy levels described in section [22] The 
format of the energ y level data file is similar to that used by 
[Harris et alJ (|2006l ). a extract from the file is given in table 

m 

In table [T] the column labelled index is the index num- 
ber given t o the c orresponding H^^CN energy level by 
[Harris et alJ (|2002a) . J and P are the exact quantum num- 
bers of angular momentum and parity, n is the number of 
the energy level in the J-P symmetry block, Eai is the value 
of the calculated energy level, iso labels the state as either 
H"CN (iso=0) or HN^^C (iso=l), ui, V2, Z, and vs are the 
approximate quantum numbers, Eiab is the lab/empirical 
energy, labell is a single character label which is either 'e' 
for a laboratory determined energy level, 'c' for an interpo- 
lated energy level or 't' for a corrected ab initio energy level, 
finally label2 is a single character label which identifies the 
calculated energy as either I for an interpolation made us- 
ing the fitted H^^CN/HN^^C rotational constants or C an ah 
mitio H^^CN/HN^^C energy from the Harris et al,: ( 2002b ) 
energy level list. The final energy level list is organised into 
symmetry blocks according to the J, P quantum numbers. 

The error on the energy level takes three forms, for a 
lab determined energy level this is the compound error of 
the line frequency measurements used to derive the energy 
level. For an energy level computed by interpolation of lab- 
oratory data this is the standard deviation on the fit of eq. 
(2] For a corrected ab initio energy level the error is the dif- 
ference between the energy predicted by the fit of eq. [5J to 
lab determined energy levels, and the corrected calculated 
H"CN or ab initio H^^CN energy. 



2.4 Partition function. 

It is essential to know the temperature dependent rotation 
vibration partition function to compute line intensities at 
thermodynamic equilibrium. The rotation-vibration parti- 
tion function of H^^CN/HN^^C was calculated by direct 
summation over all the energy levels in the energy level list, 
using: 



0-(r)=^(2J, + l)exp(-Ji) 



(3) 



where J is the angular momentum quantum number of the 
state i, Ei the energy of the state, k is the Boltzmann con- 
stant and T is temperature. For temperatures between 500 
and 10 000 K the partition function was least squares fit 
with the function: 



\og{Qrv{T)) = ^a,{logw{T)y 



(4) 



—0.93943933. The summation used laboratory determined 
energy levels where available and in preference to the calcu- 
lated H^^CN energy levels. The calculation was augmented 
with H^^CN energy levels for vibrational energies greater 
than 10 000 cm~^. At 296 K the partition function was cal- 
culated to be 152.93. This can be compared with a value 
of 148.72 calculated at 296 K for H^^CN using a similar 
procedure (jBarber et al.ll2002l ). 



2.5 Einstein A coefflcients and Laboratory 
determined band dipoles. 

The intensities of a few bands of H"CN and HN^^C 
have been measured IjSmith et al.l 1 19891 : iMaki et al.l Il995l : 
iDevi et al. 2003, 200J), these are the stretching fundamen- 
tals, the first CN stretch v^ hot band, the 2i;i, 2i;i -|- va and 
2w2 overtone bands. The laboratory determined band dipoles 
of H^^CN and H^^CN, t ogether with ah initi o band dipoles 
for H^^CN calculated bv lHarris et all (|2002al lbh are listed in 
table El The band dipoles for H"CN match those of H^^^N 
to within 15%, except for the «3 fundimental and its hot 
band. This verifies that the H^^CN Einstien A coefflcents 
are a reasonable approxomation to those of H^^CN. The vs 
bands for both H^^CN and H^^CN is unusually weak and has 
an intensity structure which shows an unusual double peak 
in the R branch. These occur at slightly different J" for 
each isotopologue. To more accurately account for the un- 
usual intensity structure of the V3 fundamental in the H^^CN 
linelist, we have u sed the band dipol es and Herman- Wallace 
constants given bv lMaki et al.l (|l995r ) to compute Einstein A 
coefficients for lines of the CN stretch fundamental and hot 
bands. Throughout this work these lab determined Einstein 
A coefficien ts are substituted for the ab initio Einstein A co- 
efficients of lHarris et al.l l|2002bl ). Einstein A coefficients for 
individual rotation vibration lines are calculated from the 
band dipole and Herman- Wallace constants with the follow- 
ing formula. 

. „ „ 3FhlFhw 2 
Aif = Cf TTTTTT-rvM 



(2J' + 1) 



(5) 



where Fhl and Fhw are the Honl-London and Herman- 
Wallace factors as described bv lMaki et al.l l|l995r) , fj, is the 
band dipole, u is line frequency, and C = 647r''/(3c^/i). To 
return an Einstein A coefficient in s~^ with dipole in Debye 
and frequency in cm~^, then C should be set to 3.136186 x 
10-^ 



3 THE NEW LINELIST 



with n = 4. The standard deviation from logjQ{Q(T)) is 
0.0075 and the coefficients are ao = —57.453223, ai — 
86.387042, 02 = -46.668113, as = 11.018178 and 04, = 



Using the Harris et al.l (12002b|) Einstein A coefficients, the 
iMaki et al.l ( 19951 ) laboratory intensity measurements, the 
lab/empirical, computed and ah initio H^'^CN/HN^'^C en- 
ergy levels, we have generated a new H^^CN/HN^'^C linelist. 
In this linelist the weak transitions between high energy 
states with vibrational energy great er than 10 000 cm~^ are 
are accounted for by using only the [Harris et al.l (|2002a ) ah 
initio data. For the majority of transitions between states 
of lower energy the calculated H^^CN/HN^^C energy levels 
are used to g ive line frequ e ncy an d the ab initio Einstein A 
coefficients of lHarris et al.l l|2002bl ) to calculate the intensity. 
However, where available, the lab/empirical H^^CN/HN^'^C 
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Table 1. A sample from the lab/empirical energy level list, which is available in full, in electronic form, 
from either the CDS archive (http://cdsweb.u-strasbg.fr/cgi-bin/qcat7/MNRAS/) or from our website 
|http://www.tampa.phys.ucl.ac.uk/ftp/astrodata[|7 
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Table 2. Available laboratory determined band dipoles (Debye) for H^'^CN, the corresponding theoretical 
and laboratory data for H^^CN are shown for comparison. Fitting errors in the last quoted digit are given 
in parentheses, where available. 



K,4.''.«3) 


{v'l,v-,V',v-) 


H"CN 


H12CN 


Hl^CN theory 


(0,2,0,0) 


(0,0,0,0) 


0.047(4) 


0.0496(2) 


0.0479(11) 


(0,0,0,1) 


(0,0,0,0) 


0.000309(2) 


0.001362(4) 




(0,1,1,1) 


(0,1,1,0) 


0.00263(3) 


0.001794(3) 




(1,0,0,0) 


(0,0,0,0) 


0.085(5) 


0.0831(17) 


0.0853(16) 


(2,0,0,0) 


(0,0,0,0) 


0.00731 


0.00881(12) 


0.0086(4) 


(2,0,0,1) 


(0,0,0,0) 


0.00065 


0.00068(1) 


0.000677(4) 



energy levels (see section 12.31 ) are used to calculate line fre- 
quency in preference to the calculated energy levels. For the 
special cases of the U3 fundamental and first hot band we 
h ave used the band dipoles and Herman- Wallace constants 
of lMaki et al.l (|l995l ) to compute Einstein A coefficients and 
intensity for individual lines. 

We have truncated the H^^CN/HN"C linelist at a 
minimum intensity of 3 x 10~^* cm molecule"^ at 3000 K. 
This results in a linelist of 34.1 million lines which accounts 
for more than 99.9% of the opacity of the full linelist at 
3000 K. Th e form at of the linelist is identical to that of 
[Harris et al.l (J2006) and a sample of the linelist is given 
in table |3] here v is frequency, E" is lower state energy, 
Ail is the Einstein A coefficient, J" and J' are lower state 
and upper state angular momentum quantum numbers, p 
is parity where 1 is even and odd, n is the number of 
the level in the J-parity block, index is the unique label 
of the energy level, iso labels a H^^CN state if or a 
HN^'^C state if 1, «i, V2, I, vz are the approximate quantum 
numbers. Where the approximate quantum numbers have 
not been assigned a value of —2 is given. The full ver- 
sion of the linelist is available from either the CDS archive 
|http://cds web.u-strasbg.fr/cgi-bin/qcat?/MNRAS/j , or from 



our website http://www.tampa.phys.ucl.ac.uk/ftp/astrodata 



IT 



Note these site provide the data in the form of table 
or alternatively as a file whose three columns are the 
wavelength in A, g„fnm.,n > m for absorption and the 
lower state energy, E" , eV. 

The untruncated linel ist can be obtained by download- 
ing the original set of i Harris et al.l (l2002b) Einstein A co- 
efficients, the new assigned energy level list and running 
the supplied FORTRAN utility program dp sort-H13CN- 
v2.1.f90 . The Einstein A coefficient file from [Harris et al.l 
l|2002bt ) is sorted using ah initio H^^CN/HN^^C frequencies, 



so that the generated H"CN/HN^^C linelist wiU no longer 
be in full frequency order. 

Figure [l] shows the intensity and frequ ency of the 
H"C N lines from the HITRAN04 database (|Rothman et all 
l2005f l against those from our new H^'^CN linelist, for the ws 
(CN stretch) fundamentals and 3112 (bend) overtone. There 
is very good agreement between line frequencies and good 
agreement between line intensities. For the 3u2 overtone our 
intensities are around 20% stronger than those from HI- 
TRAN, the agreement between intensities is better for the 
W3 fundamental. 

Figures [2 El and H show the opacity of H"CN, deter- 
mined in this work, compared to that of H^^CN computed by 
[Harris et al.l (|2006l ) at a temperature of 3000 K. The opac- 
ity has been calculated by convolving the individual lines 
with Gaussians of half-width at half-maximum of i^/4000, 
where v is the wavenumber of the line, this gives a resolu- 
tion comparable to that achieved with the SWS spectrom- 
eter on the ISO satellite. The bands shown in these figures 
all show strong Q-branches, which makes the identification 
of some bands possible at high temperature and low reso- 
lution. The identification of absorption by H^'^CN from the 
vi — V2 bands is likely to be made difficult by the strong 
background of H^'^CN. However, the V2 and V2 + V3 bands 
look to be more promising candidate for identification in a 
C-star. This is because the opacity of H^'^CN at the peak of 
the Q-branches of these bands is up to 70% stronger than 
the H'^'^CN opacity, making them more easily identifiable 
against a strong background of H^'^CN absorption. 



4 SYNTHETIC AND OBSERVED SPECTRA 
FOR WZ CAS. 

Synthetic spectra for the carbo n star WZ Gas were calcu- 
lated with the WITA6 program (|Pavlenkoll2000l ). In addition 
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Figure 1. The integrated intensity and frequency of lines of the H^^CN v^ (CN stretch) fundamental and the 3112 (be nd) overtone at 
296 K . The Unes from the linelist of this work are shown with positive intensity and Unes from the HITRAN database I IRothman et alj 
|2005^ are shown with negatively intensity. 
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Figure 2. The absorption cross-section per molecule for the H CN and H ''CN V2 + V3 (bend and CN stretch) combination bands at 
3000 K. 
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Figure 3. The absorption cross-section per molecule for the H^^CN and H^'^CN vi — V2 (bend and HC stretch) combination bands at 
3000 K. 
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Table 3. A sample from the H^^CN/HN^^C linelist, the strongest 34.1 million of which are available in electronic form, from either the CDS 
archive http://cdsweb.u-strasbg.fr/cgi-bin/qcat7/MNRAS/ or from our website http://www.tampa.phys.ucl.ac.uk Quantum numbers are: 
J rotational level, p parity, n state number in J, p block, {v\ , v^ , v'i) vibrational labels. Upper state denoted by ' and lower state by " . ind 



gives the unique index number of each state; 
lb for laboratory determined frequency. 



iso gives the isomer (0 = H '^CN, 1 = HN '^C); Ibl is ai for ab initio calculated frequency and 
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Figure 4. The absorption cross-section per molecule for the H^'^CN and H^''CN V2 (bend) fundamental and hot bands at 3000 K. 



to the same approximations, opaci ties and input par ame- 
ters as in our earlier calculations (JHarris et al.l I2OO6I ) , we 
have also included our new H CN/HN C opacity. These 
synthetic spectra are calcula ted using the best fi t model at- 
mosphere for WZ Cas from [Harris et al.l (|2006r ) , which has 
Toff= 2800 K, log(Nc/No) = -0.003, log(p) = 0.0. Syn- 
thetic spectra computed with H"CN/H^^CN ratios of 0, 
1/4 and 1 compared to th e ISO/SWS observed spectrum of 
WZ Cas ()Aoki et al.lll998l ) are shown in figure H The effect 
of H^^CN absorption can be seen in the synthetic spectra 
with H^^CN/H^^CN ratios of 1/4 and 1, however below a 
ratio of 1/4 it becomes difficult to identify H^^CN absorp- 
tion on the strong background of H^^CN lines. Although, the 
strength of the some of the H^'^CN features in the theoretical 
spectrum indicates a useful ^^C/^^C ratio maybe measured 
around 3.62 and 3.64 ^m there remain other unidentified 
opacities in this region, e.g., the strong observed bandhead 
at 3.63 ^m. 

Another interesting spectral r egion to study HCN fea- 
tures we find around 14 /^m (see [Harris et al.l l|2006l )) for 
more details. The predominant features in this region are 
the Q branches of the and ( At;2 = 1) bands. Co mputed spec- 
tra for the 2800/0.0 model atmosphere given bv lHarris et al.l 
l|2006rl are shown in Fig. [S] To simplify presentation of our 
results only a short spectral range is shown here. As we see 
from the comparison of computed spectra with and without 
a contribution from HN^'^C they differ significantly, at least 
at 14.05, 14.11, 14.23 /im, where the strong features created 
by H^^CN bands are located. This region can be considered 



as very promising for the future investigations of the carbon 
isotopic ratios in atmospheres of carbon stars. 

We have searched the ISO archive for suitable high res- 
olution data. The available data in the 14 /im regime for 
potentially suitable targets IRAS 15194-5115, RY Dra, T 
LYR, WZ CasA, Y Cvn covers a relatively large parameter 
space of stellar properties with high ^^C/^^C ratio although 
none of the available spectra offer a combination of resolu- 
tion, signal-to-noise and other features which we understood 
well enough to allow us to convincingly identify H^ CN or 
HN^'^C features. Such analysis will have to await both better 
spectra and a better characterisation of the other absorbing 
species in the same wavelength region. 



5 CONCLUSION. 

We present a new set of ab initio rotation-vibration energy 
levels for H^'^CN and HN^^C which were calculated for an- 
gular momenta of J = 0, 1, 2, 3, 5, 10, 20, 30, 40, 60 and 
for both even (e) and odd (/) parity. The states were com- 
puted up to an energy of at least 10000 -I- B{J{J + 1)) cm~^ 
above the H^^CN ground state, where B ~ 1.5 cm^^ is 
the HCN rotational constant. The quantum number assign- 
m ents were made usi ng a method similar to that described 
in [Harris et~ai] (|2006D . 

The new linelist has been incorporated into our compu- 
tations of C-rich synthetic spectra. The detailed analysis of 
the infrared spectra of C-giant star with high ^^C/^^C ratios 
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Figure 5. Synthetic and observed spectra for WZ Cas, the synthetic spectra have been displaced form the observed spectrum by 0.06 
flux units. The synthetic spectra are shown for H^^CN to H^-^CN ratios of 0, 1/4 and 1. 
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Figure 6. Synthetic spectra computed of WZ Cas computed for the pure H^-^CN line list (dashed line) and H^^CN + H^'^CN lists. In 
the last case ^'^C /^'^C = 3 was adopted. The spectra are convolved to correspond to a resolution of 10000. 



ought to take account of H^'^CN and HN^^C species. More- 
over in many cases HCN spectra probably provides the best 
chance of determining the ^^C/^^C ratios in atmospheres of 
the coolest stars as the CO bands at 2.3 /im are usually 
saturated and other molecular bands are severely blended. 
Our upgraded opacity sources can be used for the determi- 
nation of carbon isotopic ratios in atmospheres of carbon 
stars. The most promising regions are those around 3.6 and 
14 /.im where the the V2 + Ws (bend and CN stretch) com- 
bination bands and V2 (bend) fundamental and hot bands 
of H^^CN and H^^CN molecules are located, respectively. 
However, to do this requires to use the spectral data of the 
proper quality. 



Finally, the use of H^^CN and H^^CN lines for numeri- 
cal analysis of infrared spectra of evolved stars is restricted 
by the incompleteness of presently available opacity sets, 
see I Tennyson et all (|2007l ). Special attention needs be paid 
to the computation of other lines for polyatomic molecules 
such as C3, NH3, CH4 and C2H2, and their isotopologues. 
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